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Abstract 

New results on the production of charged pions in p+p interactions are presented. The data 
come from a sample of 4.8 million inelastic events obtained with the NA49 detector at the 
CERN SPS at 158 GeV/c beam momentum. Pions are identified by energy loss measure- 
ment in a large TPC tracking system which covers a major fraction of the production phase 
space. Inclusive invariant cross sections are given on a grid of nearly 300 bins per charge 
over intervals from to 2 GeV/c in transverse momentum and from to 0.85 in Feynman x. 
The results are compared to existing data in overlapping energy ranges. 



*) Corresponding author: Siegfried. Wenig@cern.ch 



1 Introduction 



The NA49 collaboration has set out to explore an extended experimental programme 
concerning non-perturbative hadronic interactions at SPS energies. This programme covers ele- 
mentary hadron-proton collisions as well as hadron-nucleus and nucleus-nucleus reactions. It 
is aimed at providing for each type of interaction large statistics data samples obtained with the 
same detector layout which combines wide acceptance coverage with complete particle identifi- 
cation. It is therefore well suited for the comparison of the different processes and to a detailed 
scrutiny of the evolution from elementary to nuclear hadronic phenomena. In the absence of 
reliable theoretical predictions in the non-perturbative sector of QCD, it is one of the main aims 
of this study to provide the basis for a model independent approach to the underlying production 
mechanisms. This approach has to rely on high quality data sets which cover both the full phase 
space and a variety of projectile and target combinations. 

The present paper addresses the inclusive production of charged pions in p+p collisions. 
In the multiparticle final states encountered at SPS energies, single particle inclusive measure- 
ments cover only the simplest hypersurface of a complex multidimensional phase space. Only 
moderate hopes may be nursed to learn enough from such measurements alone to experimen- 
tally constrain the non-calculable sector of QCD, and additional studies beyond the inclusive 
surface constitute indeed an important part of the NA49 programme. 

However, as will be shown below, the experimental situation even in the most simple case 
of inclusive cross sections is far from being satisfactory. Notwithstanding a sizeable number of 
preceding efforts especially in the SPS energy range which in part date back several decades, it 
has not been possible to obtain sufficiently precise and internally consistent data sets covering 
the whole available phase space. Precision in this context may be defined as the possibility to 
establish from the existing data an ensemble of cross sections which will be consistent within 
well defined error limits. From the physics point of view, this ensemble should permit the study 
of the evolution of inclusive yields from elementary to nuclear interactions. This evolution is 
characterized for pions by deviations of typically some tens of percent in comparison to the 
most straightforward superposition of elementary hadronic interactions. Its interpretation and 
especially any claim of connection with new physics phenomena has to rely completely on com- 
parison with elementary data. Due to the poverty of available data sets in this sector, situations 
have arisen where data from nucleus-nucleus collisions are more complete and precise than the 
elementary reference. 

It seems therefore mandatory and timely to attempt a new effort in this field with the aim 
of providing an improved elementary data base. 

This paper is arranged as follows. In Section 2, the experimental situation in the SPS 
energy range is described. Sections 3 to 6 present the data samples obtained by NA49 and a 
detailed discussion of event cuts, particle identification, cross section normalization and applied 
corrections. The final data are presented and tabulated in Sections 7 and 9. A detailed compar- 
ison to existing data including a complete statistical analysis is given in Sections 8 and 9. A 
short discussion of the results in relation to forthcoming supplementary publications terminates 
the paper in Section 10. 

2 The Experimental Situation 

A band of beam momenta extending from 100 to about 400 GeV/c corresponding to 
an interval in y/s from 14 to 27 GeV may be defined as SPS/Fermilab energy range. In this 
interval, quite a number of experiments have published inclusive particle yields [1-12]. They 
may be tentatively divided into three different categories. 



1 



The first class concerns bubble chamber applications, here mostly restricted to relatively 
small data samples of order 10 000 events each obtained with the Fermilab 30 inch bubble 
chamber. It is characterized by the absence of particle identification as far as pions are con- 
cerned. 

The second class covers spectrometer measurements with small solid angle devices mea- 
suring typically well delineated and restricted ranges in production angle and featuring complete 
particle identification. 

A third class contains large solid angle spectrometers like the EHS and OMEGA facilities 
at CERN, with the restriction that very few, if any results concerning inclusive data have been 
published from these experiments. 

In the framework of the present publication the data without particle identification and 
therefore with heavy assumptions about kaon and baryon yields for pion extraction have been 
discarded from comparison. 

We are interested in the available measurements of the double differential cross section 
of identified pions 

(i) 

dxpdp\ 

as a function of the phase space variables defined in this paper as transverse momentum and 
reduced longitudinal momentum 

PL 



Xp 



(2) 



where Pl denotes the longitudinal momentum component in the cms. 

A search for such data with pion identification practically eliminates the first and third 
class. The phase space coverage of the remaining spectrometer experiments is characterized 
in Fig. [TJl for the fixed target configurations. The very limited range of these experiments is 
immediately apparent from this plot. Data are scarce or lacking completely in the regions of 
Pt below 0.3 GeV/c and above 1 GeV/c as well as xf below 0.2. In addition it will be shown 
in Section 8.2 below that the most copious data set of Johnson et al. [6| is afflicted with large 
systematic deviations which makes it unusable for quantitative reference. 

It has therefore been decided to also include some of the rather extensive data from ISR 
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Figure 1 : Phase space coverage of existing data in the a) SPS/Fermilab and b) ISR energy range 
compared to c) NA49. Each set of symbols represents a set of measurements. 
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experiments [13-23] into the comparison in the ^fs range from 23 to 63 GeV as indicated in 
Fig. [TJx As a detailed study of s-dependence over the complete kinematic plane is outside the 
scope of the present paper, only data at x F > 0.3 have been used for the comparison. In this 
region Feynman scaling is expected to hold on a few percent level whereas a more involved 
s-dependence is present at low xp. 

The phase space coverage of the NA49 apparatus is shown for comparison in Fig.d;. This 
acceptance allows for the first time a detailed study of the low p F region and is only limited by 
statistics at p T above 2 GeV/c. In addition, there is a restricted loss of acceptance at xp > 0.7 
in the lower p T range for n + due to the interaction trigger used (see Section 3.3). 

The main aim of the present paper is to contribute new data covering the accessible phase 
space as densely and continuously as possible in a single experiment. This aim needs first of 
all a very high statistics event sample. In addition, high quality particle identification has to be 
performed over the accessible range of kinematic variables. As a third basic requirement the 
absolute normalization and the systematic error sources should be controlled on an adequate 
level with respect to the statistical uncertainties. 

3 The NA49 Experiment 

NA49 is a fixed target experiment situated in the H2 beam line at the CERN SPS accel- 
erator complex. It uses a set of large Time Projection Chambers (TPC) together with two large 
aperture Vertex Magnets (VTX1,2) for tracking and particle identification. A schematic view of 
the detector is shown in Fig. |2] with an overlay of tracks from a typical mean multiplicity p+p 
event. The details of the detector layout, construction and performance are described in [24J. 

In order to introduce and explain the data taking and analysis policies adopted in the 
present work on elementary hadronic cross sections, some general features of the detector sys- 
tem have to be stressed: 

- The data flow from the 1.8 ■ 10 5 TPC electronics channels whose analog outputs are 
digitized in 512 time buckets per channel produces event sizes of about 1.5 Mbyte after 
zero suppression. This in turn limits the total number of events that can be recorded, 
stored and processed. The data sample analyzed in this paper corresponds to a raw data 
volume of 10 Tbytes. 

- The readout system which has been optimized for operation with heavy ion collisions 
only allows a recording rate of 32 events per accelerator cycle. 

- This readout rate can be saturated with rather modest beam intensities of about 10 4 /s 
which keeps the rate of multiple events during the 50 /is open time of the TPC tracker 
small enough to be readily eliminated off-line. 

- Due to the limitations discussed above, the unbiased running on a beam trigger alone 
with alternating full and empty target is not feasible. The experiment has therefore to be 
operated with an interaction trigger which unavoidably introduces a certain trigger bias. 
The corrections for this bias are discussed in detail in Section 5.2 below. 

3.1 The Beam 

The secondary hadron beam was produced by 400 GeV/c primary protons impinging on 
a 10 cm long Be target. Secondary hadrons were selected in a beam line set at 158 GeV/c mo- 
mentum with a resolution of 0.13%. The particle composition of the beam was roughly 65% 
protons, 30% pions and 5% kaons. Protons were identified using a CEDAR ll2"5ll Ring Cerenkov 
counter (see insert Fig. El) with a misidentification probability of less than 10~ 3 . Typical beam 
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Figure 2: NA49 detector layout and tracks of a typical mean multiplicity p+p event. The open 
circles are the points registered in the TPC's, the dotted lines are the interpolation trajectories 
between the track segments and the extrapolations to the event vertex in the LH 2 target. The 
beam and trigger definition counters are presented in the inset. Detector name abbreviations are 
explained in the text. 



intensities on target were 3 • 10 4 particles per extraction of 2.37 s. Two scintillators (SI and 
S2 in Fig. |2j) provided beam definition and timing, together with a ring-shaped veto counter 
(V0) reducing the background from upstream interactions. Three two-plane proportional cham- 
bers (BPD1-3 in Fig. 0) with cathode strip readout measured the projectile trajectory yielding 
150 yum position and 4.5 /irad angle resolution at the target position where the beam profile had 
a full width at base of 6 mm horizontally and 4 mm vertically. 

3.2 The Target 

A liquid hydrogen target of 20.29 cm length (2.8% interaction length) and 3 cm diameter 
placed 88.4 cm upstream of the first TPC (VTPC-1) was used. The exact target length is de- 
termined from the distribution of the reconstructed vertex positions in high multiplicity empty 
target events as shown in Fig.|3] 

Runs with full and empty target were alternated. As the empty/full event ratio can be 
reduced from initially 18% to 9% after suitable offline cuts, and as the empty target rate has 
been shown to be treatable as a small correction to the full target data, see Section 5.1, the 
empty target running could be kept to a fraction of a few percent of the total data taking. 
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Figure 3: Distribution of reconstructed vertex position in high multiplicity empty target events 
delivering the exact target length. 



The target was filled with para-Hydrogen obtained in a closed-loop liquefaction system 
which was operated at 75 mbar overpressure with respect to atmosphere. This results, at the 
mean atmospheric pressure of 965 mbar measured over the running periods of the experiment, 
in a density of pn = 0.0707 g/cm 3 . In addition, the cross section determination has to take into 
account the density of gaseous hydrogen p ET present in the empty target runs. This density is 
obtained from the empty over full target ratio of high multiplicity events observed in a small 
fiducial volume around the target center. The density ratio Pet/ph turns out to be 0.5% which 
indicates an increase of the average temperature in the empty target by about 40°K above liquid 
temperature. The boiling rate in full-target condition has been verified to present a negligible 
modification of the liquid density. 



3.3 The Trigger 

Beam protons were selected by the coincidence C-S1-S2-V0 yielding the beam rate 
Rbeam- The interaction trigger was defined by a small scintillation counter (S4 in Fig. |3) in 
anti-coincidence with the beam. This counter of 2 cm diameter was placed on the beam trajec- 
tory 380 cm downstream of the target, between the two Vertex Magnets. The trigger condition 
C-S1-S2-V0 ■ S4 delivered the interaction rates Ret and Ret for full and empty target oper- 
ation, respectively. The corresponding trigger cross section a tr i g is calculated according to the 
formula 

P H 

at " 9 ~ ph-I-Na/A ' (3) 

where /, Na, A and p# denote, respectively, the target length, the Avogadro constant, the atomic 
number and the liquid hydrogen target density. The interaction probability in hydrogen P H 
has to be extracted from the measured rates Rft, Ret and Rbeam- Taking into account the 
exponential beam attenuation in the target, the reduction of beam intensity due to interactions 
upstream of the target, the reduction of the downstream interaction probability in full target 
operation, and the gaseous hydrogen content of the empty target cell as discussed above, it is 
determined by the relation 

pH _ Rft — Ret q Rft ~ Ret Ret Pet ^ ^ 

Rbeam 2 ■ Rbeam Rbeam PH 
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In this expression higher order terms are neglected. The resulting trigger cross section, aver- 
aged over three running periods in the years 1999, 2000 and 2002, amounts to 28.23 mb. From 
this value the total inelastic cross section as measured by NA49 is derived using a detailed 
Monte Carlo calculation which takes into account the measured inclusive distributions of pro- 
tons, kaons and pions as well as the contribution from elastic scattering in order to determine 
the loss of events due to produced particles hitting S4. The different components resulting from 
this calculation are presented in Tabled 



loss from p 
loss from it, K 
contribution from a e i 


28.23 mb 
3.98 mb 
0.33 mb 

-1.08 mb 


predicted o inei 


31.46 mb 


literature value 


31.78 mb 



Table 1 : Contributions derived by a detailed Monte Carlo calculation to the determination of the 
inelastic cross section a ine i. 

It appears that the extracted inelastic cross section compares with the literature value of 
31.78 mb ll26l to within one percent. The total trigger loss concerning inelastic events amounts 
to 14.4%. The rejected events come mostly from target diffraction (about 2.5 mb) and from 
non- diffractive events containing forward charged particles at high xf which hit S4. In order to 
account for this event rejection, a topology-dependent correction to the inclusive cross sections 
is applied, as described in Section 5.2. 

3.4 Tracking and Event Reconstruction 

The tracking system of the NA49 detector comprises a set of large volume Time Projec- 
tion Chambers covering a total volume of about 50 m 3 . Two of them (VTPC-1 and VTPC-2 in 
Fig. |3) are placed inside superconducting Vertex Magnets with a combined bending power of 
9 Tm. The magnets define by their aperture the phase space available for tracking. Two larger 
TPC's (MTPC-R and MTPC-L in Fig. El) are positioned downstream of the magnets in order to 
extend the acceptance to larger momenta and to provide sufficient track length for precise parti- 
cle identification via ionization energy loss dE/dx measurement. The use of the same detector 
for the complete range of interactions implies a separation of the sensitive TPC volumes with 
respect to the beam trajectory. Otherwise, in the case of heavy ion collisions non-interacting 
beam particles would create excessive chamber loads and the density of secondary tracks would 
be prohibitively high. The corresponding acceptance loss for low p T particles in the x F range 
above 0.5 has been remedied by the introduction of a small TPC on the beam line in between 
the two magnets (GTPC in Fig. |2) which, in combination with two strip-readout proportional 
chambers (VPC-1 and VPC-2 in Fig. 12) ensures tracking acceptance up to the kinematic limit. 
As the target is placed upstream of the magnets there is a gradual reduction of acceptance at 
low pt in the backward hemisphere which for pions starts to be effective for xp < —0.05. The 
overall acceptance region for tracking of pions in the forward Xp/pr plane has been shown in 
Fig.Ct. 

The event reconstruction proceeds through several steps. Firstly, all charged track can- 
didates leaving at least 8 space points (clusters) in the TPC system are pattern recognized and 
momentum fitted. This step includes the formation of global tracks from track segments visi- 
ble in different TPC's (see Fig. El). Secondly, a primary event vertex is fitted using all global 
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tracks found in the event together with the measured beam track. After a successful vertex re- 
construction, a second path of momentum fitting is performed also including tracks which are 
only recorded in the downstream TPC's outside the magnetic field. In this final momentum fit, 
the determined vertex point is used as an additional measured point on the track. All tracks with 
an acceptable x 2 of the fit are retained for further analysis. 

3.5 Event Selection 

Two steps of event selection are introduced in order to clean up the event sample and to 
reduce the empty target background. 

Firstly, cuts on the beam position close to the target are performed. The position of the 
incoming beam particle is registered in three Beam Position Detectors (BPD1-3 in Fig. El). 
The three corresponding measurements are required to be well defined and collinear in both 
measured transverse coordinates by imposing that the extrapolation from BPD-1 and BPD-2 
to BPD-3 coincides with the beam position measured in BPD-3, as demonstrated in Fig. HJ 
A further, less restrictive cut is imposed on the beam profile. These cuts are bias-free as the 
measurement takes place before the interaction. 




-0.4 -0.2 0.2 0.4 

Extrapolation from BPD-1 and 2 [cm] 



Figure 4: Beam position at BPD-3 versus extrapolation of beams from BPD-1 and BPD-2 to 
BPD-3. Only events falling between the two lines are accepted. 

Secondly, the interaction vertex is constrained to a fiducial region around the target po- 
sition by applying a cut on the longitudinal coordinate of the reconstructed vertex position. 
Though this cut is very efficient in reducing the empty target contribution, it must be carefully 
designed as the precision of the reconstructed vertex position depends strongly on the event 
configuration. The longitudinal vertex cut is therefore performed depending on the track mul- 
tiplicity and on the track topology. Short or very small laboratory angle tracks are not entering 
the track sample determining the vertex and the cut boundaries are placed so that no event from 
the target is rejected. This is exemplified in the normalized vertex distributions from full and 
empty target shown in Fig. and b. For the 5% of the target events with at least one recon- 
structed track, for which the reconstruction software does not give a reliable longitudinal vertex 
position, the target center is used. 

The combined event selection described above reduces the relative empty target yield 
from 18% to 9% retaining 85% of the liquid hydrogen target events. This final sample corre- 
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Figure 5: Normalized vertex distributions from full and empty target events with selected track 
multiplicity a) one and b) five and more. The imposed vertex cuts are indicated by arrows. 

sponds to 4.8 • 10 6 events which were obtained in three running periods between the years 1999 
and 2002, as shown in Table El 

3.6 Track Selection 

A track in the NA49 detector is defined by the ensemble of clusters (points) in three 
dimensions that a charged particle leaves in the effective volume of the TPC system. These 
clusters have a typical spacing in track direction of about 3 cm in the VTPC's and 4 cm in the 
MTPC regions. The tracking in a TPC environment has several decisive advantages which make 
the track selection a relatively straight-forward and safe task: 

- The local efficiency for cluster formation is practically 100% with a loss rate below the 
permille level. This is guaranteed by the choice of operation point of the readout propor- 
tional chambers and by the high reliability of the readout electronics which has a fraction 
of less than 1CT 3 of missing or malfunctioning channels. 

- The total number of points expected on each track can be reliably predicted by detector 
simulation with the exception of some edge regions for the VTPC's due to magnetic field 
E x B effects. 



Year 


Events taken 
Full target Empty target 


Events after cuts 
Full target Empty target 


1999 


1 211 k 


41.2 k 


906 k 


13.7 k 


2000 


2 648 k 


47.8 k 


2 049 k 


16.9 k 


2002 


2 508 k 


69.0 k 


1 814k 


21.8k 


Total 


6 367 k 


158.0 k 


4 769 k 


52.4 k 



Table 2: Data samples analyzed. 
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- Regions of 100% acceptance can be readily defined in each kinematic bin Axp, Ap T , 
and the azimuthal angle wedge A$ by inspecting the distribution of points per track in 
comparison to the expected value. In practice this is achieved by adjusting A$ such that 
this distribution does not show tails beyond a well-defined average. This gives at the same 
time an experimental handle to stay away from edge regions which show a drop of the 
number of points. 

- The only possibility of track losses or a reduction of track length is due to weak decays 
or hadronic interaction in the detector gas. The policy adopted in the present analysis is 
to allow a track to be shorter than expected if it has a minimum length and if the lost 
points are concentrated at the end of the track. It has been verified by eye- scans on such 
shortened tracks that either a decay (presence of one additional downstream track at an 
angle to the primary track) or a hadronic interaction (several additional tracks emerging 
from a defined interaction vertex) is present. 

In practice each track entering the analysis sample has to have at least 30 points in order to 
ensure a minimum quality of particle identification (see Section 4). This corresponds to detected 
track lengths in excess of 90 cm. The only exception to this criterion is in the extreme forward 
direction where the GTPC/VPC combination has only 9 space points and where no particle 
identification via dE/dx is feasible. An example of a point number distribution in a typical 
analysis bin is given in Fig.0 



CO 
CD 



LU 



2500 



2000 



1500 



1000 



1 p T = 0.25GeV/c 



500- 



x =0.05 

F 



LA 



50 100 150 
Number of points 



200 



Figure 6: Distribution of the number of measured points in a typical analysis bin at Xp = 0.05 
andp T = 0.25 GeV/c. 

The corresponding tracks span 3 TPC's (VTPC-1,VTPC-2,MTPC) with an expected 
number of about 160 clusters. 2.15% of all tracks are found in two distinct accumulations 
around 80 and 30 points, corresponding to tracks detected up to the end of VTPC-2 (1.7%) and 
VTPC-1 (0.45%), respectively. This has to be compared to the expected fraction from nuclear 
interactions, 1.2% and 0.3%, confirming that the reconstruction software does not introduce 
short tracks, and therefore potential background, of unexplained sources. The small fraction of 
tracks falling below the 30 point cut is corrected for by the absorption correction described in 
Section 5.4. 
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The track selection criteria defined above have been controlled by extensive eye-scans. 
This method is made very efficient in the low multiplicity environment of p+p collisions by 
the photograph-like picture built up by the measured space points and by the excellent pat- 
tern recognition capability of the human eye. Based on extensive studies with special care for 
potentially problematic areas (high p?, edge regions of acceptance, short tracks) the tracking 
efficiency is found to be 100% with an upper error limit of 0.5% below xf = 0.3 rising to less 
than 2% close to the kinematic boundary approached by the GTPC/VPC combinations. 

3.7 Acceptance Coverage, Binning and Statistical Errors 

The event sample defined above contains a total of about 28 million pions. In order to 
cover the available acceptance in an optimum fashion, a binning scheme presented in Fig.|7]has 
been chosen in the variables xp and pr- There are several aspects determining this choice: 

- Optimum exploitation of the available statistics 

- Definition of bin centers at user-friendly and consistent values of ij? and px 

- Compliance with the structure of the inclusive cross sections 

- Sufficiently small bins in p tot for optimum extraction of ionization energy loss for particle 
identification 

- Avoidance of overlaps and minimization of lost regions 

- Optimization of bin sizes for minimum binning effects and corresponding corrections 

The resulting statistical precision per bin is also indicated in Fig. |7] This precision is 
superior or equal to all other existing measurements in the SPS/Fermilab and ISR energy regions 
with the exception of a few points at large p T and/or large xp. 
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Figure 7: Binning scheme in xf and pr together with information about the statistical error. 
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4 Particle Identification 



4.1 Identification Method, Parametrizations and Performance 

The NA49 detector offers a powerful combination of tracking and particle identification 
via ionization energy loss (dE/dx) measurement in the TPC system. This system combines four 
large TPC volumes and features track lengths of between 1 m and 6 m in the kinematic region 
for pions covered in this publication. Each track is sampled by readout pads of 2.8 and 4 cm 
length in the Vertex and Main TPC's, respectively. This results in the pattern of sample numbers 
N s per track as a function of Xp and px shown in Fig. [8]. 
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Figure 8: Lines for constant number of dE/dx samples per track in the Xplpr plane (restricted 
A<£> range of selected tracks). 

The dE/dx measurement is achieved by forming a truncated mean of the 50% smallest 
charge deposits sampled along each track. The truncation transforms the Landau distribution of 
the samples into a Gaussian distribution of the mean per track if the number of samples stays 
above about 30. Under this condition the relative resolution of the energy loss measurement can 
be parametrized as 



a (N s , dE/dx) 



<7n „ (dE/dx) 



dE/dx ^ 

where the dependence of the truncated mean on the number of samples dE /dx(N s ) and the 
parameters <7 , ct and (3 are determined experimentally. To this end, for each of the 62 readout 
sectors of the TPC system, tracks are binned in laboratory momentum pi a b and a preliminary 
truncated mean distribution is formed with the samples found outside the sector under study. 
Sharp cuts are performed on this preliminary dE / dx measurement in order to separate electrons, 
pions, kaons and protons in each momentum bin. Samples from different, identified tracks in 
the given sector are then combined into truncated means with an arbitrary number of samples 
in order to determine the dependences defined above. 

The dependence of the dE/dx measurement on the number of samples, dE/dx(N s ) is 
shown in Fig.|5^. The dependences of the resolution on dE/dx and N s are presented in Fig. 
and|5£. The parameters a and (3 are fitted to 

a = -0.39 ±0.03 and (5 = 0.50 ± 0.01 . 

The NA49 TPC system employs different gases based on Ne in the VTPC and Ar in 
the MTPC. The function dE / dx(N s ) and the parameters a and (3 are found to be independent 
of the type of gas. Also the parameter a referred to the same pad length is the same for Ne 
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Figure 9: a) Percentage deviation of the mean dE / dx(N s ) from dE/dx(oo) as a function of N s ; 
resolution as a function of b) dE/dx and c) N s . All 3 panels are for MTPC. 

and Ar. This has been shown by Lehraus et al. in 1982 E71 and can be understood by the 
different behaviour of secondary ionization for different noble gases ll28ll . The only difference 
of resolution observed in the Vertex chambers is due to the shorter pad length with a dependence 
which is again given by the parameter a. The resulting values are 

^Vertex = Q41 and ^Main = a352 

The measured truncated means have a nonlinear relationship to the Bethe Bloch function 
which describes the primary ionization loss and has a different dependence on particle velocity 
((3j) in the two gases. This difference which is of the order of 0-5% for different regions of 
has to be taken into account. Using the above methodology the energy loss functions are 
independently determined for the two gases. In order to combine samples from the VTPC and 
MTPC on the same track two independent truncated means are formed. The value from the 
MTPC is transformed to the corresponding value of the VTPC using a linear transformation. 
The weighted average of these numbers, taking the respective resolutions into account, results 
in the final dE/dx measurement. 

It has to be mentioned that before the formation of truncated means a number of cor- 
rections have to be applied to each ionization sample. This concerns detailed time dependence 
extending over the complete data taking period of several years including pressure dependence 
up to second order terms, and corrections for track angles, effects of magnetic field and drift 
length dependences principally induced by the electronics threshold cut combined with electron 
diffusion in the gas. 

The resulting relative resolution is typically on the 3-4% level over most of the phase 
space covered by this experiment with tails up to 8% at low xf and pr and at large xf due 
to the decrease of sample numbers. The scatter plot of dE/dx values (referred to minimum 
ionization) versus track momentum shown in Fig. [TD| gives an impression of the performance 
with respect to the necessity to separate particles in the region of the relativistic rise. The lines 
shown in Fig. [TO] represent the mean response of the detector to the different types of particle 
derived from the measured energy loss. 

4.2 Fit Procedure and Yield Determination 

The particle identification, i.e. the determination of the yields of individual particle 
species is achieved by a fit to the dE/dx distribution for a small region of momentum defined 
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Figure 10: Energy loss dE/dx with respect to minimum ionization as a function of track mo- 
mentum p lab . 



by the analysis bin. This distribution is a superposition of Gaussians with variances following 
for each particle type the parametrization specified in Formula |5] above, and thus taking proper 
account of the variation of N s over the bin. Two typical examples are shown in Fig.[TTl 




Figure 11: dE/dx distribution for two different xf/pt bins. The line represents the fit. 



The fit has in principle nine parameters: four particle yields, four mean dE / dx values and 
the do parameter describing the width of the individual distributions. The last five parameters 
would be determined exactly by the energy loss function and the resolution function (Formula 
13) if these quantities could be absolutely predicted. Due to the complexity of the primary and 
secondary ionization processes this prediction is, however, not possible on the level of precision 
needed here. The mean energy loss which is a unique function of pi a b/ m = Pi (see Fig. 10) 
has therefore to be described by a multi-parameter approximation which must be expected to 
show deviations from the measured response. Detector related imperfections in the elaboration 
of local calibration, magnetic field effects and pulse formation introduce additional deviations 
which in general violate the f3j scaling of the elementary process. The resulting pattern of 
displacements has a smooth dependence on the kinematic variables as exemplified in Fig. fl"2l 
where the deviations are shown as a function of pr at a fixed xf for pions. 

Indeed, in order to keep the systematic error of the yield extraction well below the sta- 
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Figure 12: Deviation of the measured mean dE/dx from the predicted one as a function of px 
for pions. 



tistical error in a bin, these shifts have to be and can be determined on a sub-percent level by 
performing complete nine parameter fits. The fit procedure minimizes the x 2 over the complete 
dE/dx distribution in each bin. For the % 2 definition different prescriptions representing the 
degree to which the fit is considered reliable have been tested. Comparison of the results from 
these different prescriptions reveals consistency within the expected statistical errors, confirm- 
ing that there are no visible tails or shape distortions affecting the yield extraction. The statistical 
error of the fit is calculated from the covariance matrix of the nine parameters using the stan- 
dard x 2 definition. For the determination of the pion yield which is dominating over most of 
the measured phase space, the error turns out to be equal to the square-root of the number of 
pions in each bin with the exception of the far forward region (xf > 0.3) discussed in Section 
4.5 below. This means that the fitting method itself does not introduce any further systematic 
uncertainty to be added to the purely statistical fluctuations of the number of pions in each bin. 

As each of the 589 xf/pt bins is fitted individually without imposing external constraints, 
extensive checks of the fit outputs are performed. In fact, the results have been extracted by three 
people with two different programs and the consistency has been confirmed in about 300 mutual 
cross-check bins. 



4.3 The Region of dE / dx Crossing 

Below piab of about 3 GeV/c the energy loss functions of pions, kaons and protons ap- 
proach each other in the so-called cross-over region, see Fig.[T0l This, together with a significant 
reduction of track length in the same region, see Fig. [H prohibits independent pion identifica- 
tion. On the other hand, as pions at x F = Pt = have a pi ab of 1.3 GeV/c, this region contains 
most of the inclusive pion yield and is therefore essential to complete the data set. 

In order to make the measurement possible, a reflection technique already used in bubble 
chamber experiments (l29ll without any means of particle identification in this region is em- 
ployed. At the laboratory momentum which corresponds to xf = Pt = for pions, protons are 
at x F ~ —0.3. The proton yield can therefore be extracted from forward/backward symmetry by 
fitting protons in the symmetric bin at x F ~ +0.3. The corresponding reflection can be applied 
for kaons. 

In practice, the complete pion analysis bin is reflected with proton (or kaon) assumption. 
The dE / dx distributions for the original pion bin and the corresponding proton-reflected bin is 
shown in Fig.fPl 

The reflection technique is used for p T < 0.3 GeV/c at xp = and up to x F = 0.02 
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Figure 13: dE/dx distribution in a) the pion bin (p = 1.3 GeV/c) and b) the corresponding 
proton reflected bin (p = 48 GeV/c). 

at pr = 0.05 GeV/c. The estimated systematic error from this method is below 1% since the 
proton and kaon contributions are on the 5-10% level in this kinematic area. The consistency of 
the method has been tested in bins where both the normal extraction and the reflection method 
are usable, and its reliability was confirmed. 

4.4 Electrons at High Momenta 

Above piab of about 40 GeV/c the energy loss of pions approaches the relativistic plateau 
occupied by the electrons, rendering independent extraction difficult. This difficulty is enhanced 
by the fact that the dn ratio becomes very small, typically of the order of a few permille. In order 
to be able to control this small contribution, a Monte Carlo simulation has been constructed 
using 7T° as the main source of electrons by gamma conversion and Dalitz decay. The ir° cross 
section is obtained as the average of the measured charged pions and the combined conversion 
probability is adjusted from the measurement of the dix ratio at xp = 0.05 as shown in Fig.[T4l 
In the high xp region, the ratio is then constrained to the value predicted by the simulation and 
the small excess over the prediction obtained from the fit added to the pion yield. 




0.2 0.4 0.6 0.8 1 
p T [GeV/c] 

Figure 14: e/7r ratio as a function of p T at x F = 0.05. The line represents the Monte Carlo result. 
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4.5 Constraints from Kin and p/7r Ratios at High x p 



The region of xp above about 0.3 presents specific problems for the independent extrac- 
tion of kaons necessary to obtain bias-free pion cross sections. This is due to the fact that the 
mean ionizations of the different particle types approach each other towards saturation, and that 
the average track length steadily decreases. Furthermore, the cross sections fall steeply with Xp, 
which progressively reduces the available statistics per bin. 

The combined (K~+p)/7r~ ratio is rapidly falling with xp for all pp values. This pp av- 
eraged ratio as obtained from the NA49 data is shown in Fig. [T5h as a function of Xp. As the 
dE/dx fits become rather unreliable for values of this ratio below a few percent and at low statis- 
tics, the measurement is complemented by values from other experiments 000 also shown in 
Fig.fT5h and extrapolated smoothly to zero at xp ~ 0.85. 

For K + extraction, the situation becomes difficult already above xp ~ 0.3 due to the 
preponderant p component. This leads to problems with the independent determination of the 
kaon shift with respect to the energy loss function as described in Section 4.2 above. In fact 
the fit tends to find unphysical local minima of x 2 corresponding to large displacements of the 
K + . In this area the K + shifts have to be constrained by using the values obtained for pions and 
protons. The resulting K + /n + ratios are again compared to data from other experiments IH30D 
which agree well with the results from NA49 at xp = 0.3 presented in Fig. [T5b. This confir- 
mation of consistency is necessary, as for higher x F values the K + /n + ratio was constrained 
to that from other measurements. The rapid decrease of the n + /p ratio to values below 5% at 
xp > 0.55 imposes, together with the slightly increasing K + /n + ratio in this region, a limit on 
reliable n + extraction at xp = 0.55. 




0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 1 
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Figure 15: a) (K~+p)/7r~ ratio as a function of xp measured by |0 (open circles) and [6| (open 
squares) compared to NA49 (full triangles), b) K + /n + ratio as a function of pp measured by |0 
(open circles and triangles) and (open squares) compared to NA49 (full triangles). 



5 Evaluation of Invariant Cross Sections and Corrections 

The invariant inclusive cross section is defined as 

d 3 a 

f(x F} p T ) = E(xF,pp)-—(xp,p T ) (6) 

dp 6 

where dp 3 is the infinitesimal volume element in three dimensional momentum space. 
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This cross section is approximated by the measured quantity 



t I A 3n jpf a 3\ a trig An(x F , p T , Ap 3 ) 

f m eas{XF,PT,Ap)=E(x F ,p T ,Ap)-j ; ^ — , (7) 

where Ap 3 is the finite volume element defined by the experimental bin width, a trig the trigger 
cross section discussed in Section 3.3, N ev the number of events originating from the liquid 
hydrogen target, and An the number of identified pions from the target measured in the bin 
Ap 3 . In general and for arbitrary functional shapes of f(x F ,p F ), the measured quantity f mea s 
depends on the bin width Ap 3 via E and An. 

Several steps of normalization and correction are necessary in order to make fmeas ap- 
proach f(x F ,p T ): 

- The numbers N ev and An have to be determined from the measured full and empty target 
yields. 

- The fact that a tr i g is not equal to a ine i necessitates a trigger bias correction. 

- An has to be corrected for re-interaction of produced particles in the target volume, for 
weak decay and absorption of pions on their way through the detector, and for feed-down 
from weak decay of strange particles. 

- The finite volume element Ap 3 has to be replaced by the infinitesimal dp 3 which requires 
a binning correction. 

This list defines seven corrections which will be discussed in turn below. 



5.1 Empty Target Correction 

The particle yield An/N ev can in principle be determined from separate yield determina- 
tions in full and empty target conditions for each bin with the formula 

An\ FT ~ ET 1 / / An \ FT ( An x ET 

£ 



N e J l-e\\N e J \N t _ 

where e = R ff /R F t- 

In practice the bin-by-bin subtraction would necessitate a large enough event sample 
from empty target to comply with the statistical precision of the full target data and in addition 
would give identification problems in the regions of low cross sections where fits to the energy 
loss distributions anyhow become critical. Given the constraints in data taking of the NA49 
experiment and the sizeable reduction of the empty target rate, both evoked and described in 
Section 3, a more efficient approach is used. For the achieved empty/full target event ratio of 
9%, the empty target contribution in terms of track numbers in a given bin is in fact only about 
4-5% due to the much larger yield of empty ("zero prong") events in this sample. In addition the 
empty target events are mostly of the type p+C and p+air and produce pion yields in the forward 
(projectile) hemisphere which are very similar to p+p collisions. These points make it possible 
to extract the cross sections from full target runs alone and to treat the empty target contribution 
as a small correction. Detailed studies have shown that this correction, which is about 3-4%, 
is the same for n + and 7r~, has no measurable p T dependence, and only a slight variation with 
x F . The resulting correction, defined as the ratio of the cross section measured with the proper 
full-empty subtraction to the one extracted from full target data only, is shown in Fig. [TBI 
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Figure 16: Correction factor applied to account for the empty target contribution as a function 
of xp for the average of n + and ir~ . 



5.2 Trigger Bias Correction 

As explained in Section 3.3 the interaction trigger defined by the scintillator S4 accepts 
only 85.6% of the total inelastic cross section. This bias will reflect into the measured cross 
sections via the expression 

An 



fn 



'-'trig 



N P . 



(8) 



in a non-trivial way. It will be zero for all event topologies with no hit of S4, i.e. for the far 
forward region where - once a particle is detected there - no further particle can reach S4 
by energy-momentum conservation. On the other hand, as it has been shown that particle yields 
completely decouple between the forward and the backward hemispheres [ 301, the loss of 14.4% 
should be entirely felt in the backward region. Therefore a dependence on xp as illustrated in 
Fig. [TT] may be expected, whereas short range correlations should modify this prediction in the 
forward hemisphere. 



^ 20% 




0% 



Figure 17: Qualitative expectation of trigger bias correction as a function of xf. 

Detailed correction tables are obtained experimentally by increasing the diameter of the 
S4 counter off-line and extrapolating the observed change in cross section to surface zero. The 
method resembles the technique used in transmission experiments and relies entirely on mea- 
sured quantities. This is important to note as no event generator codes can be expected to be reli- 
able at the level of precision required here. The systematic error from this method is dominated 
by the statistical error of the evaluation, typically a factor of three smaller than the statistical 
error of the extracted data. 
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The resulting trigger bias correction over the kinematic range of the experiment as pre- 
sented in Fig.fT8k and b follows indeed the expectation for p T values above about 0.6 GeV/c. 
For smaller p T a significant micro- structure in the xp dependence appears which is a reflec- 
tion of hadronic two-body correlations driven by resonance decay (see also the discussion in 
Section 10). 




Figure 18: Trigger bias correction as a function of xf at various p T for a) n + and b) n 



5.3 Re-interaction in the Liquid Hydrogen Target 

The produced particles will travel through the target material downstream of the primary 
interaction. Given the 2.8% interaction length of the target for protons, a secondary interaction 
will occur in about 1.4% of the cases for produced hadrons. The correction for this effect is 
evaluated using the PYTHIA event generator OT1 . assuming that all daughter particles from 
these secondary interactions are reconstructed at the primary vertex. As hadronic interactions 
will also produce pions, not only make them vanish, the correction factor is > 1.0 in the high xf 
region and < 1.0 in the low xp region where pion production dominates. The x F dependence 
of the correction is shown in Fig. IT9l 



5.4 Absorption in Detector Material 

The correction for pions interacting in the downstream material of the detector is deter- 
mined using the GEANT simulation of the NA49 detector. Based on eye-scans it is assumed that 
all primary pions undergoing hadronic interactions before detection are lost. This assumption 
largely simplifies the analysis, and given the small value of the correction even in the critical 
regions of phase space, it introduces only a small systematic error. The absorption correction as 
a function of xp at two pp values is shown in Fig.|20l At low p T the xp dependence exhibits 
multiple maxima which correspond to the position of the ceramic support tubes of the TPC field 
cages O. 
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Figure 20: Correction due to absorption of produced pions in the downstream detector material. 

5.5 Pion Weak Decay 

Due to the sizeable decay length of pions in the weak channel ir — > /j, + only a small 
fraction of pions will decay on their way through the detector, ranging from a maximum of 
3.6% at the smallest detected p iab dXx F and p T = 0.0 to less than 0.5% for pi ab > 10 GeV/c. 
The bulk of these decays will cause neither a loss of the particle track nor a misidentification 
problem. The dip angle of the muon track which is most critical for its reconstruction at the 
primary vertex deviates less than 1° from the pion direction at the lowest energy, and the muon 
takes on average about 80% of the pion momentum. Muons from decays in front of the detector 
in the field-free or stray field regions are therefore reconstructed at vertex. The same applies 
for decays downstream of VTPC-2. The only sample contributing to track losses corresponds 
to decays inside the magnetic field before the primary track has left the necessary 30 points for 
analysis and where the secondary track escapes reconstruction. The tracking inefficiency for 
this sample, which amounts to 1.5% of pions at the lowest detected momentum, is determined 
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to about 20% ± 10%, by making use of the total fraction of tracks with less than about 30 
points in the low momentum range. This low inefficiency was confirmed by detailed eye-scans, 
showing that secondary interactions or kaon decays are the primary source of these short tracks, 
and not pion weak decays. The resulting correction decreases rapidly with increasing laboratory 
momentum from a maximum value of 0.36% at xf and px = as shown in Fig.l2"Tl 




P,ab l QeV/C l 



Figure 21: Correction (close to 1 /p) due to the weak decay of pions. 



As far as particle identification is concerned, muons appear shifted upwards by about one 
standard deviation from the pion position in the energy loss distribution. This small deviation is 
absorbed in the independent fit of position and width of the pion peak in each bin. 



5.6 Feed-down to Pions from Weak Decays of Strange Particles 

It may in principle be a matter of discussion whether decay pions from weak decays 
of strange particles should be subtracted or counted into the total sample. In bubble chamber 
experiments with usually a small detector size, most decays escape detection via too long decay 
lengths. For decays inside the fiducial volume, secondary vertices are readily detected and the 
corresponding decay particles eliminated from the sample. For counter experiments the situation 
is less clear. In fixed target geometry the detectors are long enough to see a sizeable fraction of 
the decay daughters and vertex fitting is usually not precise enough to eliminate all secondary 
vertices. In the early collider experiments at the ISR one can start with the assumption that 
almost all decay products are counted into the track sample. In fact none of the experiments 
compared to the present data, see Section 8, has attempted a feed-down correction. As feed- 
down pions are mostly concentrated at low px and low xf the case is saved by the fact that 
these experiments usually have no acceptance in these regions. 

In view of this situation and as the present experiment has full acceptance in the critical 
areas of phase space, a complete feed-down subtraction is performed by considering all relevant 
sources, i.e. K° s , A, S + , E~, A. The main source of systematic uncertainty is the very limited 
knowledge of the corresponding production cross sections. 

The correction is determined in three steps. First, the double-differential parent distribu- 
tions are adjusted to existing data. With this input, in the second step, a decay Monte Carlo is 
used to predict yields of daughter particles in the Xp/pr bins of the experiment. These yields 
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are, in the third step, folded with the reconstruction efficiency which is obtained from a detailed 
NA49 detector simulation using complete generated events containing the appropriate strange 
hadrons. 

The p T integrated density distributions of parent particles are shown in Fig.[22]as a func- 
tion of xp. The corresponding pp distributions are extracted for K° from averaged charged kaon 
data and for A from a combined set of bubble chamber, NA49 and ISR data fl32l . 




Figure 22: pp integrated density distribution dn/dxF of parent particles contributing to the 
feed-down correction. 



The reconstruction efficiency for the daughter pions, determined from a full GEANT 
simulation of the NA49 detector using VENUS for the input distributions, is defined as the ratio 
of reconstructed vertex tracks in a given analysis bin to the daughter tracks emitted to the same 
bin. As such, this ratio can reach values above 1.0. It has been verified that the efficiency is 
only depending weakly on the input parent-particle distribution, and can therefore be used as a 
multiplicative factor. 

The average reconstruction efficiency reaches levels of up to 50% in the low xp region. It 
is rather independent of pp but shows a strong variation with the azimuthal angle which is due 
to biases in the fitting of decay products to the primary vertex. This is exemplified in Fig. |23] 
which shows the azimuthal dependence for Lambda decay to pions with pi a b = 3.2 GeV/c in 
different bins of pp. This observation leads to the restriction of the $ window in bins which 
have otherwise full acceptance for vertex tracks (see Section 3.6). 

The resulting total feed-down correction for n + and n~ is presented in Fig.[24]as a func- 
tion of xp for various pp bins. Apparently this correction is very sizeable especially at low pp 
and extends rather far up in x F with a complex overall structure. Therefore it is no question 
that great care has to be taken in comparing experiments with undefined feed-down treatment 
in these areas. 



5.7 Binning Correction 

The extraction of invariant cross sections has by necessity to be performed in finite bins of 
phase space. The measured yield is therefore the integral of the particle density over the phase 
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Figure 23 : Reconstruction efficiency for pions with pi ab = 3.2 GeV/c from A decay as a function 
of $ for various p T values. The $ region chosen for pion extraction is indicated by full symbols, 
where $ = is in the horizontal plane. 
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Figure 24: Correction of feed-down to pions from weak decays for a) ir + and b) 7r" 



space volume element. As the density distribution as a function of the azimuthal angle is flat by 
symmetry, the problem reduces to a determination of binning corrections in x F and p T . 

Approximating the variation of particle density p in the coordinate t by an expansion in 
local derivatives, 

P (t) K p(t ) + p'(t )(t - t ) + P"{t )^^- 



the measured value corresponding to the bin center t is 

1 fto+A/2 
-A/2 



Pmeas(to) = ~T / p(t)dt » P (t ) + — P "(t )A 2 

A Jto-A/2 24 
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up to second order terms, where A is the bin width. Hence the difference between the real 
particle density at to and its measured value is proportional to the second derivative of the 
density function and to the square of the bin width. This approximation holds if the difference 
does not exceed the level of a few percent. 

For the functional forms discussed here the second derivative may be approximated from 
neighbouring data points by 



/(to) 



' A lP (t 2 ) + A 2 p(ti 
Ai + A 2 




where A x = t x — t and A 2 = t — t 2 . 

The generalization of the method to the case of double differential cross sections is 
straightforward and it can be shown that the correction can be determined independently in 
the two coordinates. 

The above consideration defines the binning correction, which can be evaluated for all 
data points if statistics permits an estimation of the second derivative. The statistical uncertainty 
induced by the correction, determined by the error of neighbouring points, is about a factor of 
10 lower than the error of the data points itself. The direct application of the correction makes 
it model or parametrization independent. 

The correction is on the 1-4% level as shown in Fig. |25]for some typical pr and xp 
dependences. In fact, one of the important aspects in the elaboration of the binning scheme, see 
Section 3.7, is to keep this correction low in the high statistics regions in order to capture the 
structures of the density function in an optimum fashion. The remaining systematic error from 
neglecting higher terms is below 0.5% as verified by Monte Carlo studies. 
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Figure 25: Correction due to the binning in a) p T and b) xp. Full circles represent the correction 
for a fixed bin of Ap T = 0.1 GeV/c and Axp = 0.05, respectively; open triangles describe the 
correction for the bin sizes actually used. 
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6 Systematic Errors 

In addition to the overall normalization uncertainty, the corrections discussed in Section 5 
introduce systematic errors which are estimated by allowing appropriate error bands for the 
underlying physics inputs. The corresponding values are given in Table |3] By summing up all 
contributions, an upper limit of 4.8% for the total systematic uncertainty can be claimed. 



Normalization 


1.5 % 


Tracking efficiency 


0.5 % 


Trigger bias 


0.5 % 


Feed-down 


0.5-1.5 % 


Detector absorption 




Pion decay n — > fi + 


0.5 % 


Re-interaction in target 




Binning 


0.3 % 


Total (upper limit) 


4.8 % 


Total (quadratic sum) 


2.0% 



Table 3: Summary of systematic errors. 



It should be noted that in the high xf region where certain assumptions have to be made 
on K + /7r + and K /n~ ratios for the purpose of particle identification, additional systematic 
uncertainties of between 1 and 3% have to be added. 

Further information on the error sources is contained in Fig. |26] which gives the distribu- 
tions of each single correction and of the resulting total correction for the 589 data points of this 
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Figure 26: Distribution of correction for a) target re-interaction, b) trigger bias, c) absorption in 
detector material, d) pion decay, e) empty target contribution, f) feed-down, g) binning, and h) 
resulting total correction. 
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experiment. It appears that the individual corrections compensate to a large extent in the overall 
result so that only few points had to be corrected by more than 10%. This might give some 
support to the estimation of the total systematic error as the sum of the squares to only 2%. This 
low value has to be regarded with some caution since the authors are the first to recognize that 
this way of minimizing possible systematic uncertainties has been in the past, one of the major 
difficulties in the quantitative comparison of hadronic cross sections. Finally, only comparison 
between different independent measurements can help to understand the problem of systematic 
error propagation. This comparison will be attempted in Sections 8 and 9 below. 

7 Results 

The following chapter presents the double differential invariant cross sections for charged 
pions resulting from the event sample and from the extraction and correction procedures dis- 
cussed above. The basic numerical information is summarized in data tables. The completeness 
and the small statistical and systematic errors of these data allow the identification of a com- 
plex structure in their dependence on the kinematic variables which cannot be parametrized by 
straight-forward arithmetic expressions. A numerical interpolation scheme has therefore been 
developed. The resulting curves are presented together with the data points in p T distributions 
at fixed Xp, in xp distributions at fixed p T , and in corresponding n + /n~ distributions. Finally, 
the cross sections are also shown in their dependence on rapidity y and transverse mass m T . 

7.1 Data Tables 

The following Tables summarize the double differential invariant cross sections as a func- 
tion of Xp and pp for the binning scheme discussed in Section 3.7 above. This scheme results 
in 281 measured values for tc + (Table H|) and 308 values for ir~ (Tabled production, the dif- 
ference for the two charges being due to the limitation of n + identification in the region above 
x F = 0.55 (see Section 4.5). 

The detailed structure of the data tables results from the attempt to cover as completely 
as possible the available phase space with bins which comply with the statistical errors of the 
data spanning five orders of magnitude in cross section. 

7.2 Interpolation Scheme 

Although the data binning has been chosen to correspond to well defined bin centers in 
x p and pt, it is desirable to also offer an interpolation scheme which produces smooth overall 
xp and pt dependences for the internal extension of the data to non-measured intermediate 
values and for comparison with other data measured at different xp and pp values. Any attempt 
at such numerical interpolation meets with the problem that the dense phase space coverage 
combined with the small statistical errors of the data reveals a micro structure both in the pp and 
Xp dependences which makes a description with arithmetic parametrizations if not impossible, 
then at least incompatible with the data quality. A multistep manual interpolation method is 
therefore selected which relies on local continuity of the cross sections both in the pp and xp 
variables and which complies with the statistical accuracy of the data points. The final result 
can be controlled by evaluating the distribution of the differences between data and interpolated 
values, divided by the statistical error of each data point. This distribution should be a Gaussian 
centered at zero and with variance unity if the interpolation scheme does not introduce a bias. 
As shown in Fig.|27]this is indeed the case for the 589 data points of this experiment. 
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Table 4: Double differential invariant cross section f(xp,pT) [mb/(GeV 2 /c 3 )] for n + produced 
in p+p interactions at 158 GeV/c. The statistical uncertainty Af is given in %. 
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Table 5: Double differential invariant cross section f(xF,PT) [mb/(GeV 2 /c 3 )] for n produced 
in p+p interactions at 158 GeV/c. The statistical uncertainty Af is given in %. 
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A/Af 

Figure 27: Histogramme of the differences A between the measured invariant cross sections 
and the corresponding interpolated values (n + and n~ combined) divided by the experimental 
uncertainty Af of the data points. 

In the Figures shown in the subsequent Sections, the interpolation for n + is extended 
beyond the region of identification, xp > 0.55, by using data from experiments at ISR energies. 
This is discussed and justified in Section 8.2. 



7.3 Dependence of Cross Section on p T and x F 

The general behaviour of the invariant cross section as a function of pp at fixed xf is 
presented in Fig.|2D In order to clearly bring out the shape evolution and to avoid overlapping 
of the interpolated curves and error bars, the values at subsequent xp values are multiplied by 
factors of typically 0.5 as given in the figure caption. 

The observed p T dependence is definitely not describable by a simple overall exponential 
parametrization at any value of xp, unless one wants to introduce local exponential slope param- 
eters for small regions of pp. In addition, distinct structures are found at low pp and xf < 0.3 
and in the pp range from 0.5 to 1 GeV/c for xp > 0.25. The structure at low pp is shown in 
more detail in the linear plots of Fig. l29l where a local maximum at pp — 0.15 GeV/c is seen to 
develop in the region 0.03 < xf < 0.2. This maximum is less pronounced for 7r~ than for ix + . 

Corresponding xp distributions at fixed values of pp are presented in Fig. [2SI Again in 
order to clarify the overall shape development, the curves for the three lowest p T values 0.050, 
0.100 and 0.150 GeV/c are multiplied by different factors specified in the figure caption. 

The low p T region is magnified in the plots of Fig. I3T1 where the measured cross sections 
are displayed without multiplication. A complicated cross-over pattern with several inflexion 
points is evident for pp values below 0.3 GeV/c which is again less developed for n~ than for 

7T+. 
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Figure 28: Invariant cross section as a function of p T at fixed x F for a) n + and b) n~ produced 
in p+p collisions at 158 GeV/c. Data and lines are multiplied successively by 0.5 for ir + and by 
0.5 up to xf = 0.35 and 0.75 for x F > 0.45 for 7r~ to allow for a better separation. 



30 



70 



p 

CM 
> 

CD 

CD 



T 



T 



a) NA49 pp-^7C + X 



x = 0.0 




: b) 



T 



T 



NA49 pp-^7C X 



0.1 0.2 0.3 0.4 0.5 0.1 0.2 0.3 0.4 0.5 

P T [GeV/c] P T [GeV/c] 

Figure 29: Invariant cross section as a function of p T at fixed x F for a) ir + and b) n~ produced 
in p+p collisions at 158 GeV/c. The behaviour in the low pr region is emphasized by using a 
linear scale. 
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Figure 30: Invariant cross section as a function of xf at fixed pr for a) n + and b) n 
produced in p+p collisions at 158 GeV/c. The steps in p T are 50 MeV/c up to p T = 
0.300 GeV/c, then 100 MeV/c up to 1.5 GeV/c and finally 200 MeV/c. Data and lines for 
Pt = 0.05, 0.1, 0.15 GeV/c are multiplied by 2.25, 1.5, 1.1, respectively, to allow for a better 
separation. 
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Figure 3 1 : Invariant cross section as a function of xf at fixed px for a) 7r + and b) 7r~ produced in 
p+p collisions at 158 GeV/c. The low xf and low pr region is magnified in order to demonstrate 
the complicated behaviour of the invariant cross section in this region. Note that data and lines 
are not multiplied by factors in contrast to Fig.l30l 
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7.4 7r + /7r Ratios 



The distinct differences in the distributions of n + and n~ over phase space result in a 
complex pattern of the corresponding n + /n~ ratios. This is demonstrated in Figs. l32l and l33l 
where the p T dependence for fixed x p and the x p dependence for fixed p T are shown in separate 
panels. 

The well-known increase of the ratios with xf up to values of about 3 at xp ~ 0.5 is 
evident (see also Fig. l46b : but for the first time with these data, well defined structures at low 
(< 0.3 GeV/c), intermediate (0.5 — 0.7 GeV/c) and high p T (> 1.5 GeV/c) as well as x F around 
0.25 become visible. 
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Figure 32: Ratio of invariant cross section for n + and n as a function of pp at fixed xp. 



7.5 Rapidity and Transverse Mass Distributions 

It is customary to view double differential invariant cross sections also on the phase space 
subsurface of y and pp although these two variables are not orthogonal in momentum space. 
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The data of Tables |4| and |5] are readily transformed into y distributions at fixed values of pt and 
shown in Fig.P31 

At fixed p T the y distribution corresponds essentially to a distorted longitudinal momen- 
tum distribution. The structures observed as a function of xp will therefore also appear at the 
corresponding y values which is indeed the case. This is visible in more detail in the magni- 
fied y distributions for px < 0.5 GeV/c shown in Fig. |35l With increasing px these structures 
happen to be less pronounced as compared to the x F distributions (Fig.EB due to the different 
projection of phase space. 
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Figure 34: Invariant cross section as a function of y at fixed p T for a) n + and b) n produced 
in p+p collisions at 158 GeV/c. The steps in p T are 50 MeV/c up to p T = 0.300 GeV/c, then 
100 MeV/c up to 1.5 GeV/c and finally 200 MeV/c. 
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Figure 35: Invariant cross section as a function of y at fixed p T for a) ir + and b) n~ produced 
in p+p collisions at 158 GeV/c. The low p T region is magnified in order to demonstrate that the 
complicated behaviour of the invariant cross section in x F is also apparent in y. 
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Another frequently used variable is the transverse mass m T = ym^ + p"^. In fact, the 
shape of the m T distributions is often claimed to be close to exponential with an inverse slope 
parameter T as predicted for particle emission from a thermal source. In Fig.|36]the invariant 
cross section is presented as a function of m T — for n + and n~ produced at y — 0.0. 



_100 



o 

cnT"" 
> 
CP 

CD 

23 
E 



0.01 



0.001 




rri - m [GeV/cf] 

T jt 



rri - m [GeV/cf] 

T 71 



Figure 36: Invariant cross section as a function of m T — m n for a) 7r + and b) n produced at 

y = 0.0 

The quality of the data allows to demonstrate that the shape of the m T distributions cannot 
be described by a single exponential and makes the determination of the local inverse slope 
T(m T — m n ) as a function of transverse mass possible. As shown in Fig. 071 there is, both for 
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Figure 37: Local slope of the m T distribution as a function of m T — m n for it + and ir . 
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7r + and for tt , a smooth and continuous variation of T with m T — around a minimum of 
T = 0.155 GeV at m T - m w ~ 0.250 GeV/c 2 extending to values in excess of T = 0.200 GeV 
both at low and and high m T — m^. 

8 Comparison to other Data 

8.1 Comparison at SPS/Fermilab Energies 

As shown in Section 2 above, there are only two preceding experiments which measured 
double differential cross sections for identified pions in overlapping ranges of */s allowing for 
direct comparison. The data of Brenner et al. from the Fermilab SAS spectrometer offer 
107 data points at beam momenta of 100 and 175 GeV/c, the data of Johnson et al. |6| com- 
prise 241 data points at the three beam momenta 100, 200 and 400 GeV/c. The corresponding 
comparisons are presented in Figs. l38l l39l and l40l where the interpolated NA49 data are used 
as reference and are shown as full lines at those values of xp and pp where direct comparison 
is feasible. As the Johnson data have not been obtained at fixed x F but with respect to the ra- 
dial scaling variable xr = E cms / E^f x , their positions have to be retransformed into xf values 
depending both on pp and y/s. Therefore in this case only x p distributions at fixed p T are given. 

Inspection of Figs . |3*S1 and l3*9l reveals that the SAS data show a good overall agreement 
with NA49, whereas very sizeable systematic deviations with a complicated pt and xf depen- 
dence appear in the Johnson data (Fig.l40l. This situation is quantified in the statistical analysis 
presented in Fig.l4~Tl 

In Fig.l4lk and d the distribution of the point-by-point statistical error of the cross section 
difference is given. These errors have most probable values around 6 and 4%, respectively, 
which are governed by the statistics of the spectrometer experiments. The long upwards tails 
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Figure 38: Comparison of invariant cross section as a function of pp at fixed xp from NA49 
(lines) with measurements from J5) at 100 (full circles) and 175 GeV/c (open triangles). 
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Figure 39: Comparison of invariant cross section as a function of xp at fixed pp from NA49 
(lines) with measurements from at 100 (full circles) and 175 GeV/c (open triangles). 
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Figure 40: Comparison of invariant cross section as a function of xp at fixed p T measured 
by NA49 (lines) with measurements from at 100 (full circles), 200 (open triangles), and 
400 GeV/c (full triangles). 
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are reflecting the decrease of cross section at high xf and/or p T . Here the NA49 errors become 
comparable or exceed the ones of the comparison data. Fig. ETb and e shows the point-by- 
point cross section difference in percent with respect to the NA49 data. If this difference is 
completely defined by statistical fluctuations, its ratio to the statistical error, shown in Fig.l4lb 
and f, should obey a Gaussian distribution with variance unity. This is indeed true for the SAS 
data as demonstrated by the Gaussian fits superimposed to the histogrammes. In fact the mean 
difference over the 107 points is less than 1% with an error of about 1%. In view of the overall 
systematic errors of 2% and 7% given for NA49 and SAS respectively, this agreement has to 
be regarded as fortuitous. More importantly however, the agreement of the data within purely 
statistical fluctuations helps to exclude systematic effects between the two measurements as 
well as s-dependences beyond the level of about 5%. 

The experimental situation is less favourable with respect to the data of Johnson et al. 
Although their statistical errors are smaller than the ones of SAS (see Fig.|4TJl) already the cross 
section differences, Fig. l41e . show a very broad and off-centered distribution. The normalized 
differences, Fig. I41F . are off by nearly 3 standard deviations whilst the variance is around 4 
units. Even taking account of a mean shift of the data by 16%, which might be in agreement with 
the normalization uncertainties mentioned in their publication, this indicates major systematic 
effects which are also clearly visible in the direct data comparison of Fig.l40l 
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Figure 41 : Statistical analysis of the difference of the measurements of [ 8 1 (upper three pannels) 
and [6| (lower three pannels) with respect to NA49: a) and d) error of the difference of the 
measurements; b) and e) difference of the measurements; c) and f) difference divided by the 
error. 
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8.2 Comparison with Forward Production at ISR Energies and Extension of Data 
Interpolation for n + 

As the identification of n + suffers from the preponderant proton component for xf > 0.6, 
it is of interest to try to make use of the large set of ISR data in the forward region in order to 
extend the NA49 results on n + to the same phase space coverage as the n~ data. The CHLM 
Collaboration, Albrow et al. lfT5l fT6l [T9l l23l has published sizeable sets of cross sections in 
the region 0.3 < xf < 0.85 and 0.3 < pr < 1.2 GeV/c. These data are more abundant 
for 7r + production where samples at yfs = 31, 45 and 53 GeV |[T9l as well as independent 
measurements at y/s = 45 GeV |23| are available. For n~ only measurements at y/s = 45 GeV 
ll23l and a set of data at one fixed angle [ 16 1 have been published. 

The ISR data can be compared directly to the NA49 results in the region of overlapping 
measurements. For n~ this is presented in Fig. 02| showing various xf distributions at fixed pp. 
A detailed statistical analysis making use of a total of 111 data points is presented in Fig. l43l 
There is good agreement over the complete range of measurements with the exception of a 
series of points obtained at yfs = 23 GeV |[T6ll which are clearly above the NA49 data by about 
25% (open triangles in Fig.l42l. 




Figure 42: Comparison of invariant cross section as a function of xp at fixed px measured by 
NA49 (lines) with measurements at a) y/s = 45 GeV El and b) yfs = 23, 31, 45, 53, 63 GeV 

EE 

For 7r + the situation is more complicated and summarized in Fig. |44j The ISR data of 
E0 obtained at y/s = 31, 45 and 53 GeV show good agreement (Fig. 03k) whereas the data 
of ll23l from the same collaboration show sizeable upward shifts of about 20% in the xf region 
where they overlap with the NA49 measurements (Fig. B4*b"). On the other hand this internal 
discrepancy vanishes for higher values of xp, also seen in Fig.l44b. It has therefore been decided 
to use the combined data of [ 1 9 1 and ll23l in the appropriate x F ranges in order to extend the 
data parametrization of NA49 for n + . The statistical analysis of the [19 1 data is presented in 

Fig. la 
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Figure 43: Statistical analysis of the difference of the measurements of ll23ll with respect to 
NA49: a) error of the difference of the measurements; b) difference of the measurements; c) 
difference divided by the error. 




Figure 44: Comparison of invariant cross section as a function of a) pr at fixed Xp published by 
ED at y/s = 31, 45, 53 GeV and b) x F at fixed p T published by E3 at = 45 GeV to NA49 
measurements represented as lines. 
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Figure 45: Statistical analysis of the difference of the measurements of lfT9ll with respect to 
NA49: a) error of the difference of the measurements; b) difference of the measurements; c) 
difference divided by the error. 



All in all, the data comparison over a range of y/s from 17 to 63 GeV demonstrates 
Feynman scaling to the level of a few percent both for n + and n~ in the range of xp > 0.4. 
A more complete study of s-dependence also in comparison to lower energies and for more 
central ranges of xp is outside the scope of this paper and will be addressed in a subsequent 
publication. 

8.3 Comparison of 7r + /7r Ratios 

A comparison of the n + /n~ ratios as a function of xf both in the Fermilab and ISR 
energy ranges is given in Fig.H^Jfor various p T values. As expected from the comparison of the 
single particle cross sections of IDE], the general agreement is satisfactory with the exception 
of the ISR data at y/s = 45 GeV in the lower range of xp. 

It is however interesting to note that also the data of Johnson et al. [6| show reasonable 
agreement as far as the ratios are concerned, which is clearly confirmed by the statistical anal- 
ysis presented in Fig. |47j This indicates that the large systematic deviations found in the cross 
sections cancel out in the ratios thus allowing a further precise cross check of the NA49 data. 



9 Integrated Data 

In addition to the double differential invariant cross sections discussed above, integrated 
quantities such as invariant and non-invariant yields and particle ratios as a function of xp or 
y, first and second moments of the p T distributions and finally the total pion multiplicity are 
of interest. Such quantities are evaluated below and compared to other available data. Integra- 
tions are generally performed numerically applying Simpsons parabolic approximation to the 
interpolated data. 
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Figure 46: Comparison of n + /n~ ratios as a function of xp for various pt values measured by 
jS) (upper four pannels), @ (middle four pannels), and |23| (lower four pannels) to the NA49 
results represented as lines. 




Figure 47: Statistical analysis of the difference of the n + /n~ ratio of [6| with respect to NA49: 
a) error of the difference of the ratio measurements; b) difference of the ratio measurements; c) 
difference divided by the error. 
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9.1 p T Integrated Distributions 



The distributions of the non-invariant and invariant yields are defined as 

dn/ dx F = tc I a ind ■ y/s/2 -J f/E- dp% (9) 

F = J f ■ d P 2 T 

dn/ dy = n / a ind ■ J f • dp% 

with f = E ■ d 3 a/dp 3 , the invariant cross section. These quantities are summarized in Table |5] 
and shown as a function of x F and y in Fig.l48l 

The statistical errors are in general below the percent level with the exception of the n~ 
data at xp > 0.55. The overall experimental uncertainties are therefore completely governed 
by the quoted systematic errors. This is in particular true for the extrapolated data for n + in the 
phase space region above x F = 0.55. 

Further integrated quantities such as n + /n~ ratios, mean transverse momentum (p T ) and 
mean transverse momentum squared (p T ) are presented in Fig.|49]again as a function of x F . The 
distinct structures and in particular the differences between n + and n~ in the first and second 
moments of the p T distributions are noteworthy and will be discussed in Section 10. 
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0.0772 0.45 


0.4940 0.19 


0.3312 0.33 


2.4 


0.2525 


0.1462 


0.55 


0.4439 0.65 


0.0792 0.65 


0.4631 0.25 


0.2964 0.50 


0.1544 0.78 


0.0275 0.78 


0.4896 0.40 


0.3232 0.74 


2.6 


0.1924 


0.1034 


0.65 


0.2046 1.00 


0.0309 1.05 


0.4587 0.60 


0.2915 1.15 


0.0542 1.46 


0.00820 1.46 


0.4865 0.82 


0.3269 1.58 


2.8 


0.1420 


0.0703 


0.75 


0.0727 2.05 


0.00955 2.10 


0.4501 1.00 


0.2804 1.75 


0.0157 2.96 


0.00205 2.96 


0.4741 1.64 


0.3068 3.05 


3.0 


0.1018 


0.0463 


0.85 


0.0200 3.25 


0.00232 3.35 


0.4040 2.60 


0.2314 3.50 


0.00355 5.93 


0.00041 5.93 


0.4345 3.44 


0.2602 6.25 


3.2 


0.0698 


0.0298 




















3.4 


0.0449 


0.0179 




















3.6 


0.0255 


0.00966 




















3.8 


0.0133 


0.00446 




















4.0 


0.00621 


0.00179 




















4.2 


0.00247 


0.000566 




















4.4 


0.000755 


0.000141 




















4.6 


0.000208 


0.0000351 




















4.8 


0.0000396 


0.0000065 



Table 6: pp integrated invariant cross section F [mb-c], density distribution dn/dxp, mean 
transverse momentum (p T ) [GeV/c], mean transverse momentum squared (p T ) [(GeV/c) 2 ] as 
a function of xp, as well as density distribution dn/dy as a function of y for n + and n~ . The 
statistical uncertainty A for each quantity is given in %. 
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Figure 48: Integrated distributions of n + and n~ produced in p+p interactions at 158 GeV/c: 
a) density distribution dn/dxp as a function of xp\ b) Integrated invariant cross section F as a 
function of x F ; c) density distribution dn/dy as a function of y. 




X F X F X F 

Figure 49: a) it + /it~ ratio, b) mean p T , and c) mean p F as a function of x F for 7r + and 7r 
produced in p+p interactions at 158 GeV/c. 
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9.2 Data Comparison 



Given the scarcity and generally incomplete phase space coverage of other experiments 
the extraction of integrated quantities is liable to suffer from sizeable systematic uncertainties. 
This becomes evident by comparing the p T integrated pion yields from to the NA49 data as 
shown in Fig.|50l Although the differential data are in good agreement with each other as shown 
in Section 8.1, the use of simple exponential or Gaussian parametrizations of the px distributions 
for extrapolation into the unmeasured regions of px 1 8 1 result in systematic deviations of the 
integrated values which reach many standard deviations with respect to the given errors. 



o 

.Q 



0.01 



0.001 




Figure 50: a) Comparison of pr integrated invariant cross section F as a function of xp for n + 
and 7T _ measured by [8| to NA49 results (represented as lines); b) Deviation of the measure- 
ments of [8 1 from the NA49 results in percent. 



The EHS experiment [10-12] which is directly comparable to NA49 in terms of phase 
space coverage and potential particle identification capability and which has accumulated a 
sizeable data sample of 470k events at a beam momentum of 400 GeV/c has only published 
Pt integrated distributions of the invariant and non-invariant cross sections. The results are 
compared to the NA49 data in Fig. |3T] As shown by the percentage deviation of the invariant 
yields (Fig.lSTb). a consistent upward shift of about 12% is evident up to xp — 0.35. This shift is 
not compatible with the precise fulfillment of xp scaling exhibited by the ISR data for xp > 0.3 
in the same energy range as discussed in Section 8.2. This indicates a normalization problem 
of the EHS data. Above xp = 0.35, the EHS results show large and evidently unphysical 
systematic effects especially for the n + yields. 

It is interesting to also compare the non-invariant density distributions dn/dxp, as pre- 
sented in Fig.|521 The increase of the percentage deviations below xp ~ 0.15 is entirely due 
to the energy dependence of the functional relation between non-invariant and invariant yields, 
see equation (9), as shown by the superimposed line in Fig. l52b . This relation predicts for s- 
independent invariant cross sections a linear increase of particle density with ^fs at Xp = 0. 
This is evidently borne out by the data. A more detailed discussion of s-dependences will be 
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Figure 51: a) Comparison of p T integrated invariant cross section Fas a function of x F for n + 
and 7r~ measured by lfl2ll to NA49 results (represented as lines); b) Deviation of the measure- 
ments of [ 12 1 from the NA49 results in percent. 
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Figure 52: a) Comparison of density distribution dn/dxF as a function of %f for n + and tt~ 
measured by lfT2l to NA49 results (represented as lines); b) Deviation of the measurements of 
[fT2l from the NA49 results in percent. 
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carried out in a separate publication. 

A comparison of pr integrated n + /n~ ratios with NA49 is possible for the data of ll8l 
and fll2|, whereas mean p T and mean p\ are only available from [ 1 1 and lfl~2l . The summary 
of these comparisons is presented in Fig.l53l 

The 7r + /n~ ratios of JH1 are in fair agreement with the NA49 results whereas the data of 
[12J deviate above xf — 0.3. The first and second moments of the pr distributions show an 
upward trend also below this xf value which complies with the expected s-dependence of these 
quantities. 
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Figure 53: Comparison as a function of xp of a) n + ratio measured by lfT2l (full circles) 
and (open symbols), b) mean px of h~ measured by [10|, and c) mean p\ for n + and %~ 
measured by [ 12 1 to NA49 results (represented as lines). 



9.3 Total Pion Multiplicity 

The integration of the dn/dxF distribution presented in Table[6]yields the following total 
pion multiplicities: 

(n n +) = 3.018 
(?v-> = 2.360 

(7T + /7T-) = 1.279 

The statistical errors on these quantities are negligible compared to the overall normalization 
uncertainty of 1.5% given in Table 

9.4 Availability of the Presented Data 

The NA49 data are available in numerical form on the Web Site ll33ll as far as the tabulated 
values are concerned. In addition and in order to give access to the complete data interpolation 
developed in the context of this publication, two large sets of 7r + and n~ momentum vectors 
(5 • 10 7 each) can also be found on this site. These files are obtained by a Monte Carlo technique 
with importance sampling without weighting. By normalizing with the respective total pion 
multiplicities and the total inelastic cross section given above, distributions of cross sections in 
arbitrary coordinate representations and with arbitrary binning may readily be obtained. This 
might be found useful for comparison with production models using Monte Carlo methods 
with finite bin sizes especially in view of the importance of binning effects demonstrated in 
Section 5.7 above. 
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10 Discussion 



The purpose of this paper is the establishment of a complete and internally consistent set 
of inclusive pion cross sections at SPS energy. As such, the data may serve as a basis of com- 
parison with other hadronic interactions, in particular with collisions involving nuclei in p+A 
and A+A reactions. There is however, another important aspect of this work which is more di- 
rectly related to the status of the theoretical understanding of soft hadronic collisions. This status 
might be described as deeply unsatisfactory due to the lack of reliability and predictive power of 
present attempts in the field of non-perturbative QCD. In the presence of sufficiently consistent 
and precise experimental information, both the multitude of ad-hoc microscopic multiparame- 
ter models and the attempts at a statistical approach to hadronic interactions fail to describe the 
detailed structure of the data even on the most primitive level of single particle inclusive cross 
sections. It is outside the scope of this publication to proceed with a detailed discussion of these 
problems which has to be left to subsequent papers. Two aspects are however worth mentioning 
here. The first concerns the important local structure visible in the data, the second addresses 
the s-dependence. 



10.1 The Importance of Resonance Effects 



Distinct local structures are visible in all inclusive distributions shown in this paper, in 
the xp or y as well as in the pp dependences, in the n + /n~ ratios as well as in the moments 
of the p T distributions. These structures are significantly different for n + and n~ production. 
Similar structures have been first observed in a more restricted experiment at the CERN PS in 
1974 [ 34- 1 and discussed in relation to resonance decay. In fact, by using the measured inclusive 
cross sections for the low-mass mesonic and baryonic resonances p and A ++ such structures are 
readily predicted if the mass distribution of these resonances is properly taken into account. This 
is illu strated in Fig . |54j where the pp dependence at fixed x p , the x p dependence at fixed pt , and 
the mean pp as a function of xf for tc + originating from p and A ++ decays are shown. Clearly, 
the salient features of the data are predicted already from this very restricted set of resonances. 
In addition, the known charge dependence of mesonic and hadronic resonance production in 
p+p collisions also implies a distinct difference between n + and n~ as observed in the data. 
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Figure 54: Monte Carlo study of a) pt distribution at fixed xp, b) xp distribution at fixed pp, 
c) mean p T as a function of xp for tt + resulting from p and A ++ decay. 
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Evidently a number of questions are raised by this argumentation. It is of course known 
since a long time ll35l that the vast majority of all hadrons come from resonance decay. There 
has however been a surprising lack of attempts to assess the real consequences of this fact 
for all aspects of the inclusive distributions. In addition to the local structures discussed above 
also the full xf and pt dependences proper can be reproduced from resonance decay if an 
appropriate yield of higher mass resonances is allowed for. The past estimates [35 1 have only 
taken account of a very limited number of mostly mesonic resonances in the low mass region. 
There is however good reason to conclude that a much larger fraction of the almost 200 listed 
baryonic and mesonic resonances is produced in p+p interactions and contributes via cascading 
decay. A measurement of inclusive resonance cross sections well beyond the lowest-lying states 
therefore becomes mandatory for any progress in understanding of both the inclusive sector and 
particle correlations. 

10.2 Energy Dependence 

The problem of the s-dependence of the invariant cross sections is touched upon in this 
paper only in the very limited forward region of xf > 0.3 and in the range of 17 < y/s < 
53 GeV. This extension in y/s is necessary due to the lack of comparison data in the more imme- 
diate neighbourhood of interaction energy. It demonstrates the s-independence of the inclusive 
cross section, to a few percent level, in this area. This finding corresponds to the well-known 
hypotheses of limiting fragmentation [ 36 1 and Feynman scaling ll371 . From the physics point of 
view these hypotheses have the dual aspect of the manifestation of fragmenting excited objects 
or of a direct trace of the partonic structure of the colliding hadrons. As such, a more detailed 
scrutiny also at lower and higher energies seems indicated. This study meets, in the lower y/s 
region, with the usual problem of the scarcity of consistent data sets. At y/s values above the 
ISR energy range there is in addition, the basic experimental problem of access to the forward 
kinematic region. Already at RHIC energies, y/s = 200 GeV, the present experimental equip- 
ment is limited to xf < 0.15 at the mean transverse momentum of the produced particles. This 
limitation becomes more and more restrictive with increasing energy at pp colliders and at the 
LHC where at best some studies around xf = and in the area of diffraction are feasible. 

The zone at xp = has, on the other hand, its proper interest due to the fact that here 
the bulk of pionic multiplicity is produced and that the invariant cross section shows a steady 
increase with y/s, the "rising rapidity plateau". These effects are not touched upon in this paper 
due to the principal problem of feed-down corrections as described in Section 5.6. As the ISR 
data for example are uncorrected for strangeness feed-down, this correction has to be performed 
before any reliable conclusions can be drawn. In a more general sense, the subsequent crossing 
of flavour thresholds and the saturation of corresponding particle yields with increasing y/s, 
from SU(2) at the PS/AGS to SU(6) at the LHC, poses considerable problems both on the 
experimental and on the interpretation level. 

11 Conclusion 

A new and extensive set of inclusive cross sections for pion production in p+p collisions 
at the CERN SPS is presented and compared in detail to existing data in the corresponding 
energy range. The new data cover the available phase space with a consistent set of double- 
differential cross sections and with systematic and statistical uncertainties well below the 5% 
level. This precision allows for the first time in the SPS energy range, the observation of a rich 
substructure in the data. This structure pervades all kinds of inclusive distributions and precludes 
any attempt at simple analytic parametrizations. It is a direct manifestation of resonance decay. 
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The importance of this building-up of the inclusive particle distributions, over the full phase 
space, from the decay of resonances is again stressed, as it has profound consequences for the 
understanding of particle production in the non-perturbative sector of QCD. The new data may 
also serve as a basis for comparison with the more complex hadronic collisions involving nuclei 
in p+A and A+A interactions. Data in these areas will be provided, with similar precision, in 
upcoming publications. 
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